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ABSTRACT 

We present new near infrared spectroscopic measurements for 29 luminous high-redshift active galactic nuclei 
(AGNs) and use the data to discuss the size and other properties of the narrow line regions (NLRs) in those 
sources. The high resolution spectra have been used to carefully model the Fe II blends and to provide reli- 
able [O III] A5007, Fe II and H/3 measurements. We find that about 2/3 of all very high luminosity sources 
show strong [O III] A5007 lines while the remaining objects show no or very weak such line. While weak 
[O III] A5007 emitters are also found among lower luminosity AGNs, we argue that the implications for very 
high luminosity objects are different. In particular, we suggest that the averaging of these two populations in 
other works gave rise to claims of a Baldwin relationship in [O III] A5007 which is not confirmed by our data. 
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We also argue that earlier proposed relations of the type ^nlr°= L,' ..., , where 7?nlr is the radius of the NLR, 
are theoretically sound yet they must break down for 7?nlr exceeding a few kpc. This suggests that the NLR 
properties in high luminosity sources are very different from those observed in nearby AGNs. In particular, 
we suggest that some sources lost their very large, dynamically unbound NLR while others are in a phase of 
violent star-forming events that produce a large quantity of high density gas in the central kpc. This gas is 
ionized and excited by the central radiation source and its spectroscopic properties may be different from those 
observed in nearby, lower luminosity NLRs. We also discuss the dependence of EW(H/5) and Fe II/H/3 on 
luminosity, black hole mass, and accretion rate for a large sample of AGNs. The strongest dependence of the 
two quantities is on the accretion rate and the Fe II/H/3 correlation is probably due to the EW(H/3) dependence 
on accretion rate. We show the most extreme values measured so far of Fe II/H/3 and address its correlation 
with EW([0 III] A5007). 

Subject headings: galaxies: active - galaxies: nuclei - galaxies: Seyfert - quasars: emission lines - galaxies: 
starburst 



L INTRODUCTION 

The narrow line regions (NLRs) of active galactic nuclei 
(AGNs) have been studied, extensively, from the ground and 
from space. This component of the nucleus is spatially re- 
solved from the ground in nearby sources and HST observa- 
tions extend the range to a redshift of about 0.5. Thus, de- 
tailed NLR mappings are now available for a large number of 
sources covering a large range of luminosity and redshift (Fal- 
cke, Wilson, & Simpson 1998; Bennert et al. 2002, hereafter 
B02, Schmitt et al. 2003). 

The spectroscopic characteristics of the NLR gas have been 
studied, extensively, over several decades and high quality 
data are now available (e.g., Veilleux & Osterbrock 1987). 
The main source of excitation of the NLR gas is photoion- 
ization by the the central continuum (see review and refer- 
ences in Netzer 1990) but shock excitation must be important 
in some parts of this region, most notably is NLRs that are as- 
sociated with jet-like radio structures (Schiano 1986; Dopita 
et al. 2002 and references therein). The gas dynamics has 
been studied too with detailed results concerning the profiles 
of various emission lines and their dependence on the level of 
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ionization, the density and the dust content of the gas (e.g., 
Veilleux 1991; Nelson & Whittle 1996; Barth et al., 2001). 
A major emphasis in recent years has been the extension of 
such works to higher luminosity sources. Some such studies 
suggest that the [O III] A5007 line width is correlated with the 
stellar velocity distribution in the bulge and thus also with the 
mass of the central black hole (Nelson, 2000; Shields et al., 
2003). 

Several recent attempts to study NLRs in large samples 
of AGNs lead to apparently conflicting results. B02 ob- 
tained narrow band HST images of seven luminous radio- 
quiet Palomar-Green (PG) quasars with z < 0.5. They argued, 
on the basis of comparison with nearby less luminous sources, 
that the NLR size (radius) scales with the [O III] A5007 and 
the H/3 line luminosities roughly as /?nlr°= L^'^- The mea- 
sured NLR sizes in their most luminous sources approached 
10 kpc (throughout this work we assume Hq = 70 km s~' 
Mpc~', fim = 0.3, and D,\ = 0.7). This dependence has been 
questioned by Schmitt et al. (2003) who studied a much 
larger sample (22 Seyfert I's and 38 Seyfert 2's), albeit with 
much lower luminosity, and found Rnlr°^ L"'^^. Croom et 
al. (2002) analyzed the specti-a of ~ 22,000 AGNs from the 
2dF quasar redshift survey (2QZ) and claimed to see a de- 
crease in the equivalent width (EW) of several narrow lines 
([O II]A3727, [Ne V]A3426, and [Ne III]A3870) with source 
luminosity. They suggested that at least part of this "Bald- 
win effect" (Baldwin 1977) is due to the increase in NLR size 
with source luminosity which leads to galactic-scale dimen- 
sions in the most luminous objects. Such NLRs are likely 
to escape the system leading to AGNs with weak or no NLR 
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emission. Testing this idea for the most intense narrow line, 
[O III] A5007, was Hmited by the low redshift, and hence rela- 
tively low luminosity in the Croom et al. (2002) ground-based 
sample. 

This paper addresses the issues of "the disappearing NLRs" 
and the Fe II/H/3 ratio in high luminosity AGNs. The work 
complements the Shemmer et al. (2004; hereafter Paper I) 
study and is based on the same data set. The paper is arranged 
as follows: §|2presents the new observations, §|3shows var- 
ious correlations involving the [O III] A5007, H/3, and Fe II 
lines and §|3discusses the implications regarding the size and 
the physics of the NLR as well as the Fe II/H/3 ratio in the 
most luminous AGNs. 

2. OBSERVATIONS 

We obtained new near infrared (IR) spectroscopic observa- 
tions for a sample of 29 high-redshift, high-luminosity AGNs 
at the Anglo- Australian Telescope (AAT) in Australia and at 
Telescopio Nazionale Galileo (TNG) in Spain. The observa- 
tions and data reduction are described in Paper I where all IR 
spectra are also shown (their Figures 1-3). 

In this work we focus on the [O III] A5007 emission line 
that, when observed, is the strongest narrow line in the spec- 
trum (all numbers given in this paper refer to the 5007A 
component only). We also show and discuss various corre- 
lations involving the broad Fe II lines. The measurement of 
the [O III] A5007 line is complicated due to the presence of 
strong, broad Fe II blends in this part of the spectrum. The ac- 
curate modeling of these blends is crucial to our study of the 
Fe II spectrum as well as the measurement of the [O III] A5007 
line. Below we give a detailed description of this process and 
illustrate the results for the case of [HB89] 1346-036. 

The near-IR spectrum of the quasar [HB89] 1346-036 
(z=2.370) was observed by Mcintosh et al. (1999). The same 
observation was later used by Yuan & Wills (2003) who re- 
measured the Mcintosh et al. (1999) spectrum and used it 
in their study of the Eddington ratio in z ^ 2 quasars. Both 
Mcintosh et al. (1999) and Yuan & Wills (2003) found an 
[O III] A5007 line in this source, which was later used in sev- 
eral of their correlations. Our superior, high S/N, better res- 
olution IR spectroscopy of the source shows broad lines due 
to Hf3 and Fe II and a weak emission feature very close to 
the putative [Olll]A5007 position (Fig. ^ top curve). We 
have used the empirical Fe II emission template of Boroson 
& Green (1992; hereafter BG92), scaled to the intensity of 
the strongest iron features in order to remove those Fe II lines. 
The broadened template is shown in the diagram and the Fe II 
subtracted spectrum is given below the original spectrum. As 
seen in the diagram, the process completely removed all trace 
of [O III] A5007 emission. We applied a similar procedure to 
the spectra of all other sources and the results listed below 
are all corrected for the Fe II blends. We suspect that other 
sources in the Mcintosh et al. (1999) and the Yuan and Wills 
(2003) samples suffer from a similar problem and hence de- 
cided not to include these objects in our analysis. We note 
that Sulentic et al. (2004) also question some of the Mcintosh 
et al. (1999) measurements. 

Table Ogives a summary of the data used in this paper. It 
includes the object name (column 1), systemic redshift (col- 
umn 2), continuum luminosity defined as XL\ at rest wave- 
length 5100A (column 3) and the basic [Olll]A5007 line 
measurements (rest-frame EW in columns 4 & 5, luminos- 
ity in columns 6 & 7, and FWHM in columns 8 & 9). Tabled 
also gives the best- fit H/3 luminosity in column 10, and the Fe 
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Fig . 1 . — An example of the [O III] A5007 line measurement process for the 
z=2.370 quasar [HB89] 1346-036 we observed at TNG in 2002. The diagram 
shows the reduced calibrated spectrum with emission lines resembling the 
[O III] A5007 doublet (top cun'e) and the Fe II subtracted spectrum {middle 
curve) where no sign of [O III] A5007 is seen. The Fe II template is also 
shown (bottom cun'e). 



ll/H/3 flux ratio^ in column 1 1 . Regarding the uncertainties on 
those numbers, some of those are discussed in Paper I and the 
others, related to the [O III] A5007 line, were obtained using 
the procedure explained in Paper I as applied to this line. 

The main result, which is apparent in Table ^ is that the 
population of high-redshift, high-luminosity quasars is di- 
vided into two distinct groups. One group (22 sources) con- 
tains objects with sti'ong [Olll] A5007 fines (EW- 10-80 A) 
and with I([0 III] A5007)/I(broad H/3) similar to the ratio ob- 
served in many low luminosity type-I AGNs. The second 
group (seven sources) shows no [O III] A5007 line within 
the observational uncertainty. To obtain the upper limits on 
EW([Olll] A5007) in those sources, we assumed a "typical" 
[O III] A5007 fine with FWHM of 1000 km s"' (about the me- 
dian in our sample, see below) and looked for the weakest 
such feature that would have been detected in our spectra af- 
ter the removal of the Fe II blends. For the best S/N spectra 
(three sources), this translates to an EW which is approxi- 
mately 0.05xEW(H/3). For the other four cases the upper 
limits correspond to 0.1 -0.2xEW(H/3). The luminosities 
corresponding to these upper limits are listed in Table[2 

The division into two distinct groups of very high luminos- 
ity AGNs is further confirmed by tfie Dietrich et al. (2002a) 
observations. These authors found that two out of the six 
luminous z ~ 3.5 quasars in their sample have prominent 
[O III] A5007 lines while the remaining four had no trace of 
this line. Combining with our new data we find that out of 35 
high luminosity quasars, 24 show strong [O III] A5007 and 1 1 
others are consistent with no such line in their spectrum. 

We have checked this finding in various ways. In particular, 
we have examined the distribution of I([Olll] A5007)/I(H/3) 
which is shown in Fig.|2] The histogram is made up of a broad 
distribution centered at about 0.3 and a group of sources with 
I([Olll] A5007)/I(H/3)< 0.1. Unfortunately, the uncertainty 
on the upper limit of I([Olll] A5007) in two of the sources 
is very large and those objects bridge the gap between the 
strong and the very weak [O III] A5007 emitters. We must 

^ This ratio matches the BG92 definition of R Fe II, i.e., the ratio between 
the EW of the Fe II blends in the A4434-A4684 band and EW{H/3) 
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Table 1 . Continuum and Emission Line Measurements. 



Quasar Name 


z" 


Log AL;^ (5 100) 


EW([0 III]) 


^|o ml 


FWHM([0 III) 


iH/3 


Fe Il/H/3 








Best Fit 


Direct 


Best Fit 


Direct 


Best Fit 


Direct 


Best Fit 








[ergs-i] 


[A] 


[A] 


[ergs-'] 


[ergs-l] 


[kms-'l 


[kms-l] 


[ergs-l] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


2QZJ001221. 1-283630 


2.339 


46.26 


< 10 




< 43.59 








44.29 


1.27 


2QZJ002830.4-281706 


2.401 


46.58 


32 


31 


44.40 


44.38 


1511 


1823 


44.63 


0.37 


UM667 


3.132 


46.28 


16 


16 


43.78 


43.79 


759 


1468 


44.44 


2.08 


LBQS 0109+0213 


2.349 


46.80 


25 


31 


44.50 


44.59 


1398 


1341 


44.93 


0.14 


[HB89] 0123+257 


2.369 


46.57 


27 


27 


44.31 


44.32 


587 


532 


44.75 


0.26 


2QZ J023805. 8-274337 


2.471 


46.57 


<7 




< 43.70 








44.70 


1.57 


SDSS J024933.42-083454.4 


2.491 


46.38 


27 


27 


44.11 


44.10 


815 


394 


44.54 


<0.1 


[HB89] 0329-385 


2.435 


46.71 


20 


24 


44.30 


44.37 


491 


450 


44.82 


0.30 


[HB89] 0504+030 


2.473 


46.32 


73 


74 


44.49 


44.51 


1065 


836 


44.38 


0.49 


SDSSJ100428.43+001825.6 


3.046 


46.44 


54 


60 


44.46 


44.51 


527 


996 


44.76 


0.59 


TON 618 


2.226 


47.31 


<3 




<44.11 








45.42 


0.65 


[HB89] 1246-057 


2.240 


47.16 


<5 




< 44.14 








45.14 


1.20 


[HB89] 1318-113 


2.306 


46.89 


14 


13 


44.31 


44.28 


1903 


2303 


44.84 


<0.1 


[HB89] 1346-036 


2.370 


46.88 


<3 




< 43.72 








45.02 


0.87 


SDSS J135445.66+002050.2 


2.531 


46.49 


<3 




< 43.22 








44.52 


0.63 


UM629 


2.460 


46.56 


35 


23 


44.40 


44.23 


1413 


929 


44.77 


1.29 


UM632 


2.517 


46.54 


15 


16 


44.03 


44.05 


196 


492 


44.85 


0.35 


UM642 


2.361 


46.29 


11 


14 


43.63 


43.74 


1696 


1859 


44.53 


0.41 


UM645 


2.257 


46.31 


23 


30 


43.98 


44.09 


525 


633 


44.67 


0.10 


SBS 1425+606 


3.202 


47.38 


23 


19 


45.03 


44.96 


1382 


659 


45.37 


0.34 


SDSS J170102.18+612301.0 


2.301 


46.34 


<7 




< 43.50 








44.50 


1.06 


SDSS J173352.22+540030.5 


3.428 


47.00 


11 


8 


44.34 


44.22 


1335 


833 


44.99 


0.39 


[HB89] 2126-158 


3.282 


47.25 


13 


9 


44.64 


44.51 


1327 


1040 


45.47 


0.49 


[HB89] 2132+014 


3.199 


45.77 


59 


52 


43.83 


43.77 


1315 


1328 


44.27 


0.54 


2QZJ221814.4-300306 


2.389 


46.54 


12 


13 


43.94 


43.96 


1791 


2343 


44.54 


0.57 


2QZ J222006.7-280324 


2.414 


47.22 


13 


16 


44.63 


44.73 


1019 


908 


45.14 


0.42 


[HB89] 2254+024 


2.083 


46.45 


15 


13 


43.94 


43.86 


612 


1390 


44.78 


1.27 


2QZJ231456.8-280102 


2.400 


46.31 


15 


15 


43.79 


43.81 


1267 


1654 


44.47 


<0.1 


2QZJ234510.3-293155 


2.382 


46.32 


20 


28 


43.96 


44.11 


887 


758 


44.79 


0.82 



Note. — The methods used to obtain 'Best Fit' and 'Direct' measurements are outlined in Paper I. 
"Systemic redshift (see Paper I). 




0.3 

F /F 
[OIII] H(3 

Fig. 2. — [O III] A5007/H/3 line ratio histogram for our sources. The dark 
region indicates sources with upper limits {arrows) on the [O III] A5007 flux. 



therefore consider two hypotheses: one is a real dichotomy 
in the [O III] A5007 line intensity and the other, a continuous 
distribution in EW([0 ill] A5007) which has a long tail at very 
small EWs. 

Very high luminosity sources that are also weak 
[O 111] A5007 emitters have been found by Yuan and Wills 



(2003) in their IR study of high redshift quasars. How- 
ever, most of those objects are broad absorption line quasars 
(BALQSOs) that are known to have a weak [O ill] A5007 
line. Regarding intermediate luminosity AGNs, BG92 find 
that 26 out of the 87 sources in their PG quasar sample show 
EW([Oiii] A5007)< lOA and 12 show EW([Oiii] A5007)< 
6A. In addition, about half of the quasars in the new Sulen- 
tic et al. (2004) high-z sample have EW([0 ill] A5007)< 6 A. 
Assuming the upper limits we obtained can be translated to 
actual EW measurements, we conclude that the EW distribu- 
tion in our sample is not very different from that of BG92. 
Combining with the information about nearby Seyfert Is we 
conclude that there are very few weak [O ill] A5007 emitters 
among low luminosity AGNs, but their fraction increases to- 
wards intermediate and high luminosity. However, Lrg mi of 
the very luminous sources in our sample is two orders of mag- 
nitude larger than observed in the most luminous BG92 ob- 
jects. This has important implications for the physics and the 
structure of the NLR in these extreme cases as discussed in 
the following sections. 

3. LUMINOSITY, EQUIVALENT- WIDTH, AND SIZE 
CORRELATIONS 

3.1. Measured and predicted NLR sizes 

B02 measured NLR sizes in seven PG quasars and com- 
pared them with sizes obtained by Falcke et al. (1998) in 
nearby Seyfert 2 galaxies. Their main finding is a strong cor- 
relation between the NLR radius (/?nlr) and the [O ill] A5007 
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Fig. 3. — The i?NLR vs. Lh/3 diagram including our luminous high-z 
quasars. Filled squares are the original B02 quasar data, and the straight 
line is the B02 RNLR~iH/3 best-fit line. Empty symbols represent high-z 
quasars with measured [O III] lines, for which Snlr was inferred from the 
B02 RNLR~i|o III] relation. One-sided error bars on the empty symbols in- 
dicate the difference between best-fit and 'direct' measurements (see text). 
Arrows represent high-z quasars, for which we have an upper limit on the 
[O III] A5007 line flux and hence a limit on Snlr from the B02 relation. Note 
the two distinct groups and the enormous predicted Rnlr at high H/3 lumi- 
nosity. 



line luminosity. Their relationship (B02 Eq. 1) scaled to the 
somewhat different cosmology adopted here, can be written 
as 



^NLR — 2.1L| 



0.52±0.06 
[O m],42 



kpc, 



(1) 



where L\ 



[O III] ,42 



= Lroni]/10 ergs . This is in perfect 
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agreement with /Jnlr°^ L,' ,„, (the uncertainty on the constant 
2.1 kpc is of order 15%). For reasons that shall become ap- 
parent later, we prefer to use the equivalent relation involving 
the H/3 luminosity (B02 Eq. 3 converted to our assumed cos- 
mology) 

/?NLR = 1.15L°HfAf°'kpc, (2) 



■"H/3,42 



"\42 , 



where Lh/3. 42 = Lh/3 / 10"*' erg s 

The recent, more detailed work of Schmitt et al. (2003) use 
a sample of 60 Seyfert 1 and Seyfert 2 galaxies and discuss 
in detail the differences between the two sub-groups, the con- 
centration of the [O III] A5007 emission, etc. The main find- 
ing which is relevant to our work is the following relationship 
(adjusted to the cosmology used here) 



R 



maj 



^•^^[Om],42 "^P^' 



(3) 



where R^^^ is the size of the semi-major axis of the 
[O III] A5007 nebulosity. This is significantly different from 
the B02 results in both the R-L dependence and the normal- 
ization. Schmitt et al. (2003) have also investigated the corre- 
lation when the seven B02 sources are added to their sample. 
The result is a steeper dependence of the form 7?maj "^ ^?o ml ■ 
A new work by Bennert et al. (2004) argues that much of the 
difference is due to orientation since most of the sources in 
Schmitt et al. (2003) are Seyfert 2s while the more luminous 
sources in B02 are all type-I AGNs. We shall return to this 
issue in §|4] 
Theoretical suggestions that the NLR size should scale with 

1/2 

^ion ' where Lion is the ionizing source luminosity, have been 
discussed in many papers (see Netzer 1990 for references 




Fig. 4. — Two dimensional spectrum around the [OlII] A5007 region of 
2QZ J222006. 7— 280324. The coordinates of the vertical and horizontal axes 
are given in pixels, where each pixel in the spatial (vertical axis) corresponds 
to 0.446" (~ 3.5 kpc). Note that there is no trace of extended [O III] A5007 
emission beyond a ~ 1" (~ 7.2 kpc) radius from the center 



prior to 1990 and Dopita et al 2002 for more recent publi- 
cations). This is based on the assumption that both the broad 
line region (BLR) gas and the NLR gas are photoionized by 
a central source whose spectral energy distribution changes 
only slightly with source luminosity. Spectroscopic studies 
show a remarkable similarity between the emission line spec- 
trum of high and low luminosity AGNs. This suggests that the 
ionization parameter, U (defined here as the ratio of the Ly- 
man continuum photon density to the hydrogen number den- 
sity A^h)^ and the typical gas density, are basically the same in 

all sources. Since U ^ Lioa/NuR^, we find R °^ L^^^^ . An ad- 
ditional assumption is that most of the [O III] A5007 emission 
originates in radiation-bounded clouds, because of the fairly 
uniform value of I([Olll] A5007)/I(narrow H/3). This allows 
to replace Lion by the luminosity of any hydrogen recombina- 
tion line, e.g., H/3. We note that there is a well established 
R-L correlation for the BLR gas (e.g., Kaspi et al. 2000) 
where reverberation mappings show that Rblr °^ L^-^^^-^ , 

A real physical explanation of any R-L dependence is still 
lacking since the mechanism controlling the gas density and 
location is unknown. Some papers assume a stratified, radia- 
tion bounded NLR with a pre-chosen run of density and hence 
level of ionization (e.g., Netzer 1990; Komossa & Schulz 
1997). Such models are naturally normalized in incident flux 
units (Lioa/R^) and can be tuned to produce the same mean 
U for all sources. However, some narrow emission lines are 
probably produced in density bounded gas (e.g., Binette, Wil- 
son, & Store hi-Bergmann 1996) which considerably compli- 
cate the models. It is also clear (e.g., Alexander et al. 1999) 
that well studied NLRs contain a large range of conditions 
with a spread in ionization parameter and gas density. An 
almost orthogonal approach is provided by the "locally opti- 
mally emitting cloud" (LOC) model (Ferguson et al. 1997). 
The assumption in this case is of a large range of conditions 
(density and covering factor) at each location, where the in- 
tensities of the various lines reflect the line production effi- 
ciency at each location. This efficiency is the highest for 
densities that are close to the critical density of the line in 
question. The model provides a natural scaling of Lion with 
^NLR provided there is a large reservoir of gas with simi- 
lar properties in all AGNs and on all scales. Finally, there 
are equally complex NLR models where shock excited gas 
contributes significantly to the observed NLR emission (e.g., 
Contini, Prieto, & Viegas 1998; Schiano 1986; Dopita and 
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Sutherland 1995) but no natural R-L scaling. 

The recent papers by Dopita et al. (2002) and Groves et 
al. (2004) provide a more solid foundation to NLR modeling. 
These authors assumed dusty, stratified NLR clouds where the 
external radiation pressure acts mostly on the dust particles 
and forces the local ionization parameter to certain specific 
values such that in the [O III] A5007 producing gas, U ~ 10"^. 
The model suggests a natural ^nlr°= i'' ^ dependence. It also 
implies that the composition and temperature of the gas are 
rather different from those assumed in other models because 
of metal depletion. 

3.2. New [O III] and R^i^^ measurements 
The combination of the B02 observations and the recent 

1/2 

theoretical developments point to a "natural" Rnlr°^ I^io^ ^^' 
pendence yet raises a severe problem regarding the NLR size 
in high luminosity AGNs. Any such scaling will lead, at large 
enough Ljon, to sizes that are larger than galactic sizes. The 
B02 results addressed here are used for normalizing this rela- 
tionship, but the problem exists at high luminosity whether or 
not their scaling is correct. Given this, we would not expect 
to see any strong [O III] A5007 emitters in high luminosity 
sources, yet our new observations clearly show such objects. 

To define the problem in a more quantitative way we plot 
in Fig. |3] three quantities vs. Lh/3 (luminosity of the entire 
emission line). The first is the measured /?nlr from the B02 
sample (seven sources) where we also show the best (modi- 
fied) B02 fit (Eq.|2}. The second is from our newly observed 
H/3 and [O III] A5007 lines with two additional sources from 
Dietrich et al. (2002a). For these we use Eq. [0 to guess 
i?NLR, given the observed Liq m] . The predicted ^nlr for 
most sources in this group lie close to the value predicted 
from Lh/3 (the straight line) confirming the small scatter in 
I([Olll] A5007)/I(H/3). The third group includes the seven 
sources from our sample and the four sources from Dietrich 
et al. (2002a) where no [O III] A5007 has been detected. For 
these we use /?nlr derived from the upper limits on Liq nil ■ 
The upper limits on /?nlr obtained in this way are a factor of 
2-3 smaller than those derived from Lh/3- 

The implications of Fig.|3]are clear For those sources with 
measured H/3 and [O III] A5007 lines, the derived /?nlr is 
enormous, exceeding 70 kpc in the most luminous sources. 
We consider those sizes completely unreasonable for reasons 
that are discussed in the next section. Like Fig.|2l this diagram 
suggest a dichotomy in the properties of the high luminosity 
quasars, where some sources show strong [O III] A5007 lines 
and others show no or very weak such emission. 

There are two other ways to verify, experimentally, the 
B02 claim. Most of our sample is at z ~ 2.5. At this red- 
shift, and the chosen cosmology, the angular diameter dis- 
tance is a weak function of redshift and corresponds to about 
7 kpc per arc-sec. The predicted /?nlr> using the B02 rela- 
tionships and our measured [O III] A5007 luminosities, cor- 
responds to a total extent of 2"- 10". This can be tested by 
spatially resolved space and ground-based observations. The 
2D specti-a of two of our sources ([HB89] 0329-385, and 
2QZ J222006.7-280324) show prominent [O III] A5007 emis- 
sion which allow such measurements (see Fig. |4}. In these 
cases, most of the line flux (> 99%) is emitted within the cen- 
tral four pixels corresponding to a radius of 7.2 kpc at each 
source. This size is a strong upper limit since much of the 
flux is likely to be due to the PSF (corresponding to ~ 1" at 
the time of observations). The two upper limits on /?nlr ob- 
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Fig. 5. — Baldwin relationship for [OlII] A5007. Mean (over luminosity 
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dotted line. 



tained in this way are a factor of ^ 5 smaller than the radii 
derived from the B02 relationships. 

The term "NLR radius" used by B02 is ambiguous since 
those authors used very low surface brightness features to de- 
fine the dimension of the [O III] A5007 nebulosity. We used 
the HST archive to extract and re-analyze the B02 images. 
In particular, we examined the source showing the largest 
[Olll] A5007 nebulosity, PG 1049-005, and remeasured its 
observed [Olll] A5007 image. We found that 95% of the Une 
emission is encircled within a radius of 1.1" which corre- 
sponds to a radius of 5.5 kpc at the source. This is half of the 
radius deduced by B02 and suggests that the bulk of the NLR 
emission is emitted within a volume which is much smaller 
than inferred by their relationships. The different way of 
measuring the [O III] A5007 nebulosity is probably the main 
source of discrepancy in normalization (i.e., the NLR radius 
at Liq jjji 42 = 1 ) between the B02 and the Schmitt et al. (2003) 
works. 

In summary, our new observations contradict the B02 re- 
sults in two ways. First, about one third of the high luminos- 
ity sources show no trace of an NLR. Second, there are direct 
indications in three cases, and sound theoretical reasons (see 
below), to suggest that most of the NLR emission is restricted 
to a volume which is much smaller than inferred from the B02 
relationships. 

3.3. The Baldwin relationship for [O III] 

The new data show the presence of two groups of lumi- 
nous AGNs, those with strong [O III] A5007 and those with 
no (or very weak) such line. Out of the 35 sources investi- 
gated by us (6 from Dietrich et al. 2002a and 29 from our 
IR sample), 24 belong to the first group and 11 to the second. 
As mentioned above, earlier studies, like the Yuan & Wills 
(2003) work, already found very weak [O III] A5007 in sev- 
eral AGNs. Most or all of these are known BALQSOs while 
only two of the 10 sources discussed here (one from Dietrich 
et al. 2002a, [HB89] 0105-265 and one, [HB89] 1246-057, 
from our sample) show BALQSO properties. Thus, weak or 
no [O III] A5007 seems to be a common property of many lu- 
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Fig. 6.— FWHM([Oin]A5007) vs. (a) AL;^(5100), (b) Mm- and (c) L/Lem for the Shields et al. (2003) sample (empty symbols) and the high-z quasars 
presented in this paper (filled symbols). The strongest correlation is between FWHM([0 III] A5007) and luminosity. The correlation weakens as the dependence 
on luminosity drops from L through Mbh (°^ iP^) to L/L^n (« L"^), as indicated by the Spearman rank congelation coefficients at the top left of each panel. 



minous AGNs. 

To further illustrate this point we plot in Fig. |5] 
EW([0 III] A5007) versus ALa(5100A) for the high-z quasars 
from our sample and that of Sulentic et al. (2004), for which 
B magnitudes were transformed to ALa(5100) assuming Li, oc 
i/~" with a = 0.5. We also included the BG92 sample in the 
diagram, and the 2QZ data of Croom et al. (2002), for which 
b] magnitudes were transformed to ALa(5100) using the same 
methods as above. The diagram shows that for the popula- 
tion of strong [O III] A5007 emitters, there is no reduction of 
EW([0 III] A5007) with source luminosity. On the other hand, 
there are many weak, or no [O III] A5007 emitters at high lu- 
minosity that could give the impression that the line EW de- 
creases with increasing source luminosity. (In fact, the Croom 
et al. 2002 data represent mean EW([0 III] A5007) in several 
luminosity bins. Assuming a fraction of weak [Olll] A5007 
emitters in those bins similar to the one found by BG92, we 
find that the plotted average may underestimate the typical 
EW([Olll] A5007) in those sources by about 25%). We sus- 
pect that the Dietrich et al. (2002b) claim of a Baldwin rela- 
tionship for this line is the result of their averaging together 
four very weak [O III] A5007 emitters with two sources show- 
ing "typical" EW([0 III] A5007). 

For completeness, we tested also the correlation of 
EW([Olll] A5007) with accretion rate (in terms of the 
LBoi/i'Edd ratio, hereafter L/Ledd; see Paper I for more details) 
Correlations with accretio n rate are found to be extremely im- 
portant in Paper I and in §|331below. For EW([Olll] A5007), 
this correlation is not significant (see Table|2}. 

3.4. FWHM correlations 

We have tested our sample for possible correlations of 
luminosity, black hole (BH) mass, and accretion-rate with 
FWHM([Olll] A5007). Such coiTelations have been in- 
vestigated in the past (e.g.. Shields et al. 2003), and 
FWHM([Olll] A5007) has been suggested as a potentially 
useful surrogate for cr*. 

Fig.|6lshows FWHM([Olll] A5007) vs. source luminosity, 
Mbh, and L/Ledd for our sample and the Shields et al. (2003) 
sources. To obtain the intrinsic line width we assumed 



where AAj, 



AA^b 
is the instrumental resolution. 



AAtrue — AAobs- AAjnst, 



(4) 

Since the slit- 



width and the seeing disk sizes were comparable during the 
time of observations, this is a reasonable assumption for 
line profiles that are indistinguishable from a Gaussian. The 
uncertainty on FWHM([Olll] A5007) is quite lai-ge for the 
poorer S/N spectra and is of the order of the instrumental (or 
rebinned instrumental) resolution {^ 600 km s~'). BH masses 
and L/L^dd for all sources on these diagrams were calculated 
as prescribed in Paper I. 

On their own, the measured FWHM([0 III] A5007) for our 
high-z sources show no correlation with luminosity, BH mass, 
or accretion-rate, due to the narrow luminosity range of our 
sample. The additional Shields et al. (2003) data increase 
this range considerably and show significant correlations with 
all three parameters. The strongest correlation is between 
FWHM([Olll] A5007) and luminosity. The strength of the 
correlation increases with the dependence on luminosity, i.e. 
less significant correlations as one goes from L to Mbh (°^ 
LO'')toL/LEdd(-i°-^). 

3.5. Correlations involving Hf] and Fe II 

In Paper I and in Tabled we give L^p, which we trans- 
formed to EW(H/3), and Fe II/H/3. Here, again, we have 
tested the correlations of these quantities against L, Mbh, and 
L/LEdd using our sample, and the samples of Sulentic et al. 
(2004) and BG92 (accretion rates for the BG92 sample were 
calculated as prescribed in Paper I). For EW(H/3) we find sig- 
nificant correlations with both L and L/L^dd- Details of the 
correlations are given in Tableland in Fig.0 In the case of 
Fe II/H/3 the only significant correlation is with L/LEdd- De- 
tails of those correlations are also given in Table |2l Fig. |8] 
shows the strongest correlation of Fe II/H/3 (against the accre- 
tion rate) for our new sources combined with those of BG92. 

The Baldwin relationship found here for EW(H/3), and 
shown in Fig.0is interesting since it is in contradiction with 
the Croom et al. (2002) finding for this line. However, 
our combined sample is far from being complete, in partic- 
ular at the high luminosity range where we specifically chose 
to observe high luminosity sources. Thus it is likely that 
lower luminosity sources at high redshift would have a larger 
EW(H/3) that will spoil the correlation. Regarding the de- 
pendence on accretion rate, similar selection effects may be 
operating but the observed correlation is so strong that we 
suggest that this may be a real effect. We also note that in 
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Table 2. Spearman Rank Correlation Coefficients Matrix 



Property vs. 


Mbh 


L/LEdd 


Fe Il/H/3 


EW([0 III]) 


EW(H/3) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




XLx(5W0) 


0.85" 


0.48" 


-0.04'' 


-o.og-i 


-0.37S 




Mbh 




0.03" 


-0.22*^ 


0.14= 


-els'" 


B 


L/Lem 






0.48= 


-0.13= 


-0.47'' 




Fe Il/H/3 








-0.39f 


-0.33'' 


01 


EW([0 III]) 










0.28'' 


£ 


EW(H/3) 















Note. — Significant correlations with chance probabilities smaller than 
1% are given in boldface. Upper limits were not included in the correlations. 
"118 sources from BG92 and this work 

124 sources from BG92, Sulentic et al. (2004), and this work 
= 107 sources from BG92 and this work 

''121 sources from BG92, Sulentic et al. (2004), and this work 
= 104 sources from BG92 and this work 

'no sources from BG92, Sulentic et al. (2004), and this work 
^135 sources from BG92, Sulentic et al. (2004), and this work 
''118 sources from BG92 and this work 
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Fig. 7.— EW(H/3) vs. (a) luminosity and (h) L/LEdd for the BG92 sam- 
ple (including broad-line AGNs, empty circles, and T<lhSls filled circles), the 
Sulentic et al. (2004) sample {empty diamonds in panel a) and the high-z 
quasars presented in this paper (filled squares). 



Fig. [7] the narrow-line Seyfert 1 galaxies (NLSls) are situated 
very close to the high luminosity, high accretion rate sources. 
The situation resembles the Baskin & Laor (2004) finding that 
EW(C IV A 1549) depends strongly on L/Ledd even in low lu- 
minosity AGNs that do not show the Baldwin relationship for 
this line. 

Fig- m is hard to interpret. Considering only sources with 
measurable Fe II/H/3 (i.e., a ratio larger than ~ 0.1) we find 
a clear trend of increased Fe II/H/3 for higher accretion rates 
(see Table|2}. In particular, we note the similar location on the 
diagram of our high-z quasars and the NLS 1 s. However, there 
is a significant number of sources with very small Fe II/H/3 
and very large L/Le(J(J- We note, in this respect, that our new 
observations extend the measurement of Fe II/H/3 to values of 
L/Lgdd never investigated before. The only exceptions, per- 
haps, are a few BALQSOs in the Yuan & Wills (2003) sample 
with spectral properties very different from those of our (non- 
BAL) sources (note that these authors measured the entire op- 
tical Fe II blends and their values must be scaled down by a 
factor 3.58 for comparison with the measurements presented 
in our work; B. Wills private communication). 

The above two findings suggest that the main reason for the 
increase of Fe II/H/3 with the accretion rate is the decrease of 
EW(H/3) with L/LEdd- In fact, looking at EW(Fe II) against 
L/Ledd (not shown here) we find no correlation at all. The 
fractional increase in Fe II/H/J with accretion rate is also con- 
sistent with the decrease in EW(H/3). However, we cannot 
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Fig. 8. — Fe Il/H/3 vs. L/LEdd for the BG92 sample, and the high-z quasars 
presented in this paper (symbols are as in Fig.0. Arrows represent no de- 
tection of iron emission (in the BG92 sample the upper limit on the Fe II/H/3 
ratio was set to 0.09, which is the lowest value they reported, and in our 
sample the upper limits are given in TablelT]. 



rule out the possibility that changes in the Fe/H abundance ra- 
tio are involved too, as shown, in Paper I, to be the case for 
N/C. 

In Fig. Owe plot the well known anti-correlation between 
EW([0 III] A5007) and Fe II/H/3 (e.g., BG92) for our sources, 
for the low-z BG92 sources, and for the new intermediate-z 
sample of Sulentic et al. (2004). Our sources are consistent 
with this trend. Fe II/H/3 for our sources is generally higher 
than in BG92 and Sulentic et al. (2004) indicating, perhaps, 
the strong dependence of this property on the total luminosity. 
However, there are clear exceptions, i.e., sources with very 
high luminosity, yet relatively small Fe II/H/3. Here, again, we 
note similar large values of Fe II/H/3 in several of the BALQ- 
SOs of Yuan & Wills (2003) . 

4. DISCUSSION 
4. 1 . Do enormous NLRs really exist? 

A major goal of the present investigation is to test the 
NLR properties and the NLR spectrum in very high luminos- 
ity sources where the theoretically predicted /?nlr^ as well 
as the empirical B02 relationship, results in unreasonably 
large dimensions. The B02 sample akeady included claims 
for sources with /?nlr"^ 10 kpc and our bright [Olll] A5007 
emitters would continue this relationship to enormously large 
NLRs (more than 100 kpc in diameter, see Fig. |3}- There 
are, however, fundamental problems in this suggested size- 
luminosity interpretation on both observational and theoreti- 
cal grounds. 

B02 suggested two estimates of the NLR size (Eq. [2 & 
13. We separate the discussion of those claims into two, ac- 
cording to the two sub-groups discovered here (the weak and 
the strong [O ill] A5007 emitters). For the weak [O ill] A5007 
emitters we find a clear contradiction between the /?nlr based 
on the observed H/3 luminosity and the one based on the 
observed (or the upper limit on) [O III] A5007 luminosity 
(Fig. |3}. As for the strong [O III] A5007 emitters, 2D spec- 
tra of two of our sources rule out the large predicted dimen- 
sions (§1^). Moreover, our new measurements of the largest 
[O 111] A5007 nebula in the B02 sample are also in conflict 
with the B02 relationships. Based on the evidence in hand 
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we suspect that the /?nlr°= i^io,, dependence breaks down at 
some intermediate luminosity scale and that the "true" NLR 
radius, defined here as the radius encompassing 95% of the 
line emission, does not exceed a few kpc even in the most 
luminous quasars. 
There are other predictions that make us question the ex- 

1/2 

istence of such enormous NLRs. The suggested ^nlr°= ^jon 
relationship would predict NLR sizes, in the most luminous 
AGNs, that exceed the size of the largest known bulges and, 
in fact, the size of the largest known galaxies (except for some 
cD galaxies). Such sizes are unacceptable for several reasons. 
First, the escape velocity from a spherical galaxy is roughly 

290Mj( ^7okpc'™^ ^~'' where Mil is the mass in units of lO'^ 
Mq and Rio the radial distance in units of 10 kpc. The new 
FWHMs listed in Table [2 and shown in Fig. |6l compared 
with the predicted /^nlr. suggest therefore, dynamically un- 
bounded NLRs. Given those velocities and dimensions, the 
dynamical time for the most luminous [O III] A5007 emitters, 
is a few x 10^/?io years suggesting a short lived phenomenon. 
The inferred amount of ionized gas in such NLRs, given radi- 
ation bounded gas and a "typical" ionization parameter, is 



Mnlr ^ 10' 



0.1 



RioNii Me 



(5) 



where Cf is the covering fraction and N21 the column den- 
sity in units of 10^' cm~^. Simple photoionization argu- 
ments show that for strong [Olll] A5007 emitters N21 > 1. 
This would give Mnlr > lO' M^-, for all of our sources and 
Mnlr > 10 Mq for the most luminous [O III] A5007 emit- 
ters. Judging by the observed FWHM, most of this material 
is probably unbound and thus flows from the center at out- 
flow rates approaching 10^ R\q Mq yr~^ assuming a density 
of ~ 10^ cm~^. All those numbers seem incompatible with 
long term mass ejection in AGN-hosting galaxies. 

The Schmitt et al. (2003) results alleviate some of these 
difficulties because of the smaller size normalization and the 
flatter R^lr-Liq ml dependence. However, it is not at all clear 
that extrapolatmg their results to much higher luminosity with 
the suggested L^'^ dependence (Eq.|3 is justified in view of 
the B02 measurements. 

We are facing a situation where sound physical arguments 
(the narrow line spectrum, the Groves et al. 2004 dusty NLR 
model, etc.) support the expected Rnlr°^ L^'^ relationship, 
yet its application to the most luminous quasars give unrea- 
sonably large sizes, masses and mass outflow rates. As ex- 
plained below, the most likely explanation in our opinion is 
that "typical" NLRs (i.e., those similar in their properties to 
the ones observed in Seyfert 1 galaxies) cannot last very long 
in high and perhaps also in intermediate luminosity AGNs. 
This means that the [O III] A5007 emitting regions in our sam- 
ple may be of a different origin and physical properties. 

4.2. Luminous NLRs as star fanning regions 

If compact NLRs are indeed typical of many high luminos- 
ity AGNs, then their properties must be very different from the 
properties of those NLRs observed in nearby sources. In par- 
ticular, the gas density in the highest luminosity NLRs must 
be several orders of magnitude larger Consider for example a 
maximum NLR size of ~ 3 kpc and assume a similar ioniza- 
tion parameter in all NLRs. This would mean that the gas den- 
sity in the most luminous [Olll] A5007 emitters is 10^- lO"* 
times larger than the density in nearby Seyfert 1 galaxies. 
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Fig. 9.— EW([OlII] A5007) vs. Fe Il/H/3 for the BG92 sample, the Su- 
lentic et al. (2004) sample, and the high-z quasars presented in this paper 
(symbols are as in Fig. Q. Upper limits on Fe Il/H/3 (arrows pointing left) 
are similar to those shown in Fig.|8] Upper limits on EW([0 III] A5007), i.e., 
no [O III] A5007 detection, for the BG92 (Sulentic et al. 2004) sample were 
set to 1 ( 1 .3) A, respectively, which is the lowest value each of them reported. 
The upper limits on EW([0 III] A5007) for our sample are given in TablefTI 



The spectroscopic properties must be very different too which 
can, in principle, be tested by accurate observations. Unfor- 
tunately, present day IR spectroscopy is very limited in this 
respect because of the restricted wavelength bands available 
to ground-based observations. 

A possible origin of a high density gas in kpc-scale nuclear 
regions is violent star-forming activity. Such events can pro- 
duce high density, large column density, non-solar composi- 
tion dusty gas. The overall spectrum of such regions is likely 
to differ from the spectrum of nearby, lower density NLRs. 
At present we can only observe [O III] A5007 which is the 
strongest emission line under a variety of conditions. Future, 
space-born spectroscopy, will be able to test this idea by look- 
ing for other emission lines. 

The scenario we propose to explain the observations of 
our high luminosity AGNs, and the apparent break down of 

1/2 

the ./?NLR°^ ^ion relationship, is of two distinct populations. 
One where such scaling continues to high luminosity due 
to radiation pressure force or other effects. This results in 
short-lived enormous NLRs that will show basically no nu- 
clear narrow emission lines during most of their life. The 
other group is those sources where starburst or another un- 
known process ejects high density gas into their nuclear re- 
gion. This gas is ionized and excited by the central radia- 
tion source and produces the observed strong [O III] A5007 
lines. Such "star-forming NLRs" would have spectral proper- 
ties that are rather different from those observed in nearby 
less luminous sources. Given the similar fraction of weak 
[O III] A5007 emitters in BG92 and in our new sample, we 
suggest that the phenomenon is of a continuous nature and 
starts at some intermediate luminosity. Its clearest manifes- 
tation is in the highest luminosity sources such as the ones 
observed here. 

4.3. Correlations involving Fe II, H/?, and [O III] 

The sample used here, which combines data from various 
different sources, is not complete. However, it allows the 
best test, so far, of the BG92 Fe II/H/3 -[O III] A5007 rela- 
tionship and several other suggested correlations at the high 
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end of the AGN luminosity function. The results presented 
here suggest that the most extreme values of Fe II/H/3 require 
very high luminosity. They also suggest that higher accretion 
rate results in larger Fe II/H/9 (Fig.|8|i and that L/Ledd is the 
most important factor for determining several other correla- 
tions. There are, however, some exceptions, i.e., sources with 
large L/Lgdd yet small Fe II/H/9. Our work shows that the 
decrease of EW(H/3) with increasing L/Ledd in the present 
sample is probably the cause for the increase of Fe II/H/? with 
the accretion rate. EW(Fe II) by itself does not depend on 
the accretion rate and it is therefore possible that the iron-to- 
hydrogen abundance ratio plays no role in this correlation. On 
the other hand, we do not have the data or a good enough 
theory (e.g., Verner et al. 2004 and references therein) to 
completely rule out the possibility that the iron abundance de- 
pends on luminosity or on the accretion rate. In a similar way, 
EW([0 III] A5007) is not correlated with L/Z-Edd, but we sug- 
gest that the Fe II/H/3 - EW([0 III] A5007) anti-correlation is 
driven by the accretion rate. 

Sulentic et al. (2004) and others argued that the general 
AGN population should be divided into two sub-groups ac- 
cording to the so-called Eigenvector 1 (El, see BG92) which 
contains four main observables. According to this scheme, 
group A sources show the strongest El properties. Most 
radio-loud AGNs belong to group B, i.e., those with weaker 
El properties. Here we are not interested in El as such since 
our sample is not complete and so is the combination of our 
sample with the BG92 and the Sulentic et al. (2004) samples. 
We are however in a position to look at the extreme end of 
the distribution in L and in L/L^di and test for those quanti- 
ties that depend on them. We have akeady discussed all the 
relevant correlations. Here we note that all those involving 
L/LEdd show that NLSls occupy the same part of parame- 
ter space as the high luminosity, high accretion rate AGNs 
observed by us. These sources are thus the high luminos- 
ity analogs of NLSls and should, perhaps, be referred to as 
narrow-line type 1 quasars (NLQls). 

5. CONCLUSIONS 

We have discussed the near-IR spectra of 29 newly ob- 
served high-luminosity high-redshift AGNs and used the data 
to argue that previous claims for expected and observed 

1/2 

^NLR°= L,Q ml dependence are in conflict with the observa- 
tions. About 2/3 of all very high luminosity sources show 
strong [O III] A5007 lines while the remaining objects show 
no or very weak such line. We suggest that the NLR proper- 
ties in high luminosity, strong [O III] A5007 emitters are very 
different from those observed in nearby AGNs, and possibly 
imply denser than "usual" NLRs. The origin of the high den- 



sity gas is likely to be a violent star-forming event in the nu- 
cleus. We also investigated Fe II/H/3 and EW([Olll] A5007) 
at the high end of the luminosity and L/L^dd in AGNs and 
showed that the first of those is probably driven by the over- 
all accretion rate while the second is independent of source 
luminosity (i.e., no Baldwin relationship) or accretion rate. 
The Fe II/H/3 ratio may be an iron abundance indicator but 
this cannot be proven observationally because of the EW(H/3) 
dependence on accretion rate. 
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